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ABSTRACT 

We study the nonlinear growth of cosmic structure in different dark energy models, 
using large volume N-body simulations. We consider a range of quintessence models 
which feature both rapidly and slowly varying dark energy equations of state, and 
compare the growth of structure to that in a universe with a cosmological constant. 
We use a four parameter equation of state for the dark energy which accurately re- 
produces the quintessence dynamics over a wide range of redshifts. The adoption of 
a quintessence model changes the expansion history of the universe, the form of the 
linear theory power spectrum and can alter key observables, such as the horizon scale 
and the distance to last scattering. We incorporate these effects into our simulations 
in stages to isolate the impact of each on the growth of structure. The difference in 
structure formation can be explained to first order by the difference in growth factor 
at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. 
We find that quintessence models that are different from ACDM both today and at 
high redshifts {z ~ 1000) and which feature late (z < 2), rapid transitions in the 
equation of state, can have identical baryonic acoustic oscillation (BAO) peak posi- 
tions to those in ACDM. We find that these models have higher abundances of dark 
matter haloes at z > compared to ACDM and so measurements of the mass function 
should allow us to distinguish these quintessence models from a cosmological constant. 
However, we find that a second class of quintessence models, whose equation of state 
makes an early (z > 2) rapid transition to w — —1, cannot be distinguished from 
ACDM using measurements of the mass function or the BAO, even if these models 
have non-negligible amounts of dark energy at early times. 

Key words: Methods: N-body simulations - Cosmology: theory - large-scale structure 
of the Universe 



1 INTRODUCTION 

Determining whether or not the dark energy responsible for 
the accelerating expansion of the Universe evolves with time 
remains a key goal of physical cosmology. This will tell us 
if the dark energy is indeed a cosmological constant or has 
a dynamical form as in quintessence models. The nature 
of the dark energy determines the expansion history of the 
Universe and hence the rate at which cosmological pertur- 
bations grow. In this paper we investigate the influence of 
quintessence dark energy on the nonlinear stages of struc- 
ture formation using a suite of N-body simulations. 

The simplest candidate for dark energy is the cosmo - 
logical constant, A (see e.g. the review bv lCarrolll lioOlh. 



Despite the success of ACDM (cold dark matter and cosmo- 
logical constant model) at fitting mu ch of the available ob- 
servational data (jSanchez et al.|[2009l ). this model fails to ad- 
dress two important issues, the fine tuning problem and the 
coincidence problem. The fine-tuning problem arises from 
the vast discrepancy between the vacuum energy level pre- 
dicted by particle physics, generically given by A'', where 
A is the physics scale considered, and the value of missing 
energy density inferred cosmologically, p ~ 10~'*''GeV*. In 
the standard model of particle physics, A could be at the 
Planck scale, A ~ lO^^GeV. The coincidence problem refers 
to the fact that we happen to live around the time at which 
dark energy has emerged as the dominant component of the 
Universe, and has a comparable energy density to matter. 



* E-mail: elise.jennings@durham.ac.uk 



Quintessence models were devised to solve the fine tun- 
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ing and coincidence problems of ACDM. In these mod- 
els, the cosmological constant is replaced by an extremely 
light scalar field whi ch evolves slo wly ( Ratra & Peebles 
'l988 ; .Wetteric h 1988; Caldwell et aLflQQS : Fcrreira fc Joyce 
19981 ). An abundance of quintessence models has been pro- 



posed in the literature which can resolve the coincidence 
problem and explain the observationally inferred amount of 
dark energy. Models of quintessence dark energy can have 
very different potentials, V{(p), but can share common fea- 
tures. The potentials provide the correct magnitude of the 
energy density and are able to drive the accelerated expan- 
sion seen today. The form of the scalar field potential deter- 
mines the trajectory of the equation of state, w(z) = P/ p, as 
it evolves in time. Hence, different quintessence dark energy 
models have different dark energy densities as a function 
of time, Q.o-b{z). This implies a different growth history for 
dark matter perturbations from that expected in ACDM. 

Cosmological N-body simulations are the theorist's tool 
of choice for modelling the final stages of perturbation 
collapse. The overwhelming majority of simulations have 
used the concordance ACDM cosmology. Here we sim- 
ulate different dark energy models and study their ob- 
servational signatures. A small number of papers have 
used N-body s i mulations to test sca lar field cosmologies 
llMa et al.ll 19991: iLinder fc Jenkins|[20"0 3: Klypin ct al. 2003; 



Francis et al.l 
20091 : lAhmi et al. 



2008; Grossi fc S pringell 120091 : ICasarini et all 
2009 ). Rather than explicitly solving 
for different potentials, it is standard practice to mod- 
ify the Friedmann equation using a form for the dark en- 
ergy equation of state, w{z). Previous work used a vari- 
ety of parametrizations for w{z), the most common be- 
ing the two parameter equation, w = wo + (1 — a)wa 
l|Chevallier fc PolarskilboOll : iLinded [20031') or the empi rical 
three parameter equation proposed by IWetterichI l|2004 ) for 
the so-called early dark energy models. The disadvantage 
of using a 1 or 2 variable parametrization for w is that it 
cannot accurately reproduce the dynamics of a quintessence 
model over a wide range of redshifts. If we wish to reproduce 
the equation of state of the original scalar field to within 
5%, a two parameter equation of state will not be able to 
achie ve this precision for a wide range of quintessence poten- 
tials (|Bassett et al.l [20041 ). Instead, we take advantage of a 
parametrization for w{z) which can describe a wide range of 
different models. In this work we use a four parameter dark 
energy equation of state which can accurately reproduce the 
original w{z) for a variety of dark energy models to bet- 
ter t han 5% for redshifts z < 10^ ijCorasaniti fc Copelandl 
I2OO3I ). 

In this paper we present three stages of N-body simu- 
lations of structure formation in quintessence models. Each 
stage progressively relaxes the assumptions made and brings 
us closer to a full physical model. In the first stage, the 
initial conditions for each quintessence cosmology are gen- 
erated using a ACDM linear theory power spectrum and 
the background cosmological parameters are the best fit val- 
ues assuming a ACDM cosmology. The only departure from 
ACDM in this first stage is the dark energy equation of state 
and its impact on the expansion rat e. In the second stage , we 
use a modified version of CAMB ([Lewis fc Bridlell2002l ) to 
generate a consistent linear theory power spectrum for each 
quintessence model. The linear theory power spectrum can 
differ from the power spectrum in ACDM due to the pres- 



ence of non-negligible amounts of dark energy during the 
early stages of the matter dominated era. This power spec- 
trum is then used to generate the initial conditions for the 
N-body simulation which is run again for each dark energy 
model. The third and final stage in our analysis is to find 
the values for the cosmological parameters, Qmh^, flhh^ and 
Ho (the matter density, baryon density and Hubble param- 
eter) such that each model satis fies cosmological distance 
constraints. Recently [Xiimi et al.l (|2009, ) used CMB and SN 
data to constrain the parameters in the quintessence po- 
tential and the value of the matter density, fim/i'^, for two 
models. In this paper we allow three parameters to vary 
when fitting each quintessence model to the available data. 
This distinction is important as changes in these parame- 
ters may produce compensating effects which result in the 
quintessence model looking like ACDM. For example, for 
a given dark energy equation of state, a lower value of the 
matter density may not result in large changes in the Hubble 
parameter if the value of Ho is increased. In going through 
each of these stages we build up a comprehensive picture 
of the quintessence models and their effect on the nonlinear 
growth of structure. 

This paper is organised as follows. In Section 2 we dis- 
cuss quintessence models and the parametrization we use for 
the dark energy equation of state. We also outline the ex- 
pected impact of different dark energy models on structure 
formation. In Section 3 we give the details of our N-body 
simulations. The main power spectrum results are presented 
in Section r4.3l Intermediate results are presented in Sections 
14.11 and 14.21 which the reader may wish to omit on a first 
pass. In Section lT4l we present the mass function predictions. 
In Section 14.51 we discuss the appearance of the baryonic 
acoustic oscillations in the matter power spectrum. Finally, 
in Section 6 we present our conclusions. 



2 QUINTESSENCE MODELS OF DARK 
ENERGY 

Here we briefiy review some general features of quintessence 
models; mor e detailed descr ip tions c an be fou n d, fo r 



ample, in iRatra fc Peebles! ll 19881); IWetterich[ d 19881) 



iFerreira fc Jovcd ([l998l ): ICopeland et all (|2006l ) and lLinded 
(|2008l ). The main components of quintessence models are ra- 
diation, pressureless matter and a quintessence scalar field, 
denoted by ip. This dynamical scalar field is a slowly evolving 
component with negative pressure. This multifiuid system 
can be described by the following action 



"2k 



' df,<pd„<fi ~ V{(p)) , (!) 



where R is the Ricci scalar, £m+r is the Lagrangian den- 
sity of matter and radiation, k = SttG, g is the deter- 
minant of a spatially fiat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) metric tensor g^i/ and V{^p) is the scalar 
field potential. We assume that any couplings to other fields 
are negligible so that the scalar field interacts with other 
matter only through gravity. Minimising the action with re- 
spect to the scalar field leads to its equation of motion 



<fi + 3Hip + 



dip 



0, 



(2) 
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where H is the Hubble parameter and we have assumed the 
field is spatially homogeneous, ip{x,t) — if(t). The impact 
of the background on the dynamics of </? is contained in 
the 3-ff<^ term. The Hubble parameter for dynamical dark 
energy in a flat universe is given by 

H\z) 



^0 



(an (1 + zf + (1 - an)e^'''« + + , (3) 



where Hq and 0,-^ = Pm/Pcrit are the values of the Hubble 
parameter and dimensionless matter density, respectively, at 
redshift 2 = and Pcit = 3//o / (SttG) is the critical density. 
The dark energy equation of state is expressed as the ratio 
of the dark energy pressure to its energy density, denoted 
as w = P/p. Once a standard kinetic term is assumed in 
the quintessence model, it is the choice of potential which 
determines w as 

In general in these theories if the contribution from the ki- 
netic {ip = 0) and gradient energy {dup/dx — 0) is negligible, 
then the effect of the scalar field is equivalent to a cosmolog- 
ical constant which behaves as a perfect fiuid, with P = — p 
or w = —1. 

2.1 Classes of quintessence models 

Two broad classes of quintessence models can be used 
to solve both the fine-tuning and coincidence problems. 
The first is based on the idea of so called 'tracker 
fields' (jSteinhardt et all 1 19991 ). These fields adapt their be- 
haviour to the evolution of the scale factor and hence 
track the background d ensity. The other class is referred 
to as 'scaling so lutions' l|Halliwelllll987l : IWands et al.lll993l : 
IWetterich|[l995l ). In these models the ratio of energy densi- 
ties, Pip/pB, is constant. 

In tracking models, the ip field rolls down its potential, 
V{ifi), to an attractor-like solution. The great advantage of 
these models is that this solution is insensitive to the ini- 
tial conditions of the scalar field produced after inflation. 
A general feature of these tracking solutions is that as the 
scalar field is tracking behind the dominant matter compo- 
nent in the universe, its equation of state, w^p, depends on 
the background component as 

— fj3 (w-n-w^) 

PB 

where pB and u;b denote the background energy density and 
equation of state respectively, with wb = 1/3 (radiation era) 
and Wb = (matter era). As a result, the energy density of 
the scalar field remains sub-dominant during the radiation 
and matter dominated epochs, although it decreases at a 
slower rate than the background density. The quintessence 
field, p^, naturally emerges as the dominant component to- 
day and its equation of state is driven towards w — ~1. An 
example of a tracking m o del is the inverse potential form 
proposed bv lZlatev et all (|l999l ). V(ifi) ~ M-^+^if-", where 
M is a free parameter that is generally fixed by the require- 
ment that the dark energy density today Que ~ 0.7 and so 
the quintessence potential must be V ~ Pcrit. This implies 
that ifi is of the order of the Planck mass today, (p ~ Mpi. 
With a ^ 6, the quintessence field equation of state is ap- 
proximately Wo ^ —0.4 today. 



(5) 



In scaling quintessence models, the ratio of energy den- 
sities, p^/ pB, is kept constant, unlike tracking models, where 
p^p changes more slowly than pB. During the evolution of the 
energy density in a 'scaling' model, if the dominant matter 
component advances as p oc a~", then the scalar field will 
obey Q,ip = n'' jo? after some initial transient behaviour. 
Scaling quintessence models can suffer from an inability to 
produce late time acceleration, whilst at the same time ad- 
hering to observational constraints, such as, f or example, 
the l o wer limit on dur i ng nu cleosynthesis l|Bean et al.l 
iooj). lAlbrecht fc SkordisI (|200G| ) used a modified coefli- 
cient in their scaling potential, V{}p) = Vpe""^*^, where 
Vp((y3) = (<p — -B)" -I- A, resulting in a model which can pro- 
duce late time acceleration as well as satisfyi ng cosmological 
bound s, for a variety of constants A and B. iBarreiro et al.l 
considered a linear combination of exponential terms 
in the scalar fleld potential and found this yielded a larger 
range of acceptable i nitial energy den s ities f or compared 
with inverse models. ICopeland et~aLl l|2000l ) also consider 
supergravity (SUGRA) corrections to quintessence models, 
where the resulting potential can exhibit either 'tracking' 
or 'scaling' behaviour depending on which path the scalar 
fleld takes down its potential towards the minimum where 
it would appear as a cosmological constant. 

The physical origin of the quintessence fleld should 
be addressed by models motivated by high energy particle 
physics. As the vacuum expectation value of the scalar fleld 
today is of the order of the Planck mass, any candidates 
for quintessence which arise in supersymmetric (SUSY) 
gauge theories may receive supergravity corrections which 
will alter the fleld's potential. It is this fact that motivates 
many authors to argue that any quintessence model inspired 
by particle physics p otentials must be based on SUGRA. 
iBrax fc Martini l|l999l ) discuss such models and employ the 
potential ^(vs) = A''+"/<p"e'"/2'^'' with a value of q > 11 in 
order to drive too close to —1 today. 

In summary, in this paper we will consider six 
quintessence models which cover the behaviours discussed 
above. In particular, INVl and INV2, which are plotted in 
Fig. [TJ have inverse power law potentials and exhibit track- 
ing solu tions. The INVl model i s the 'INV model consid- 
ered by ICorasaniti fc Copelandl l|2003l ) and has a value of 
Wo = —0.4 today. As current observatio nal data favour a 
value of TOO < -0.8 jSanchez et al.ll2009l 'l. the INVl model 
will be used as an illustrative model. We shall consider a 
second inverse power law model (INV2) which i s in better 
agreement with the constraints on w. As noted bv lCorasanitil 
(2004), the scale A in the inverse power law potential, 
V'(ip) = h"^* /ip°' is flxed by the value of Q.-de today. Solving 
the coincidence problem requires this scale for A to be con- 
sistent with particle physics models. For values of a ^ 6 it is 
possible to have energy scales of A ~ 10® GeV. Setting a = 6 
results in an equation of state with wo = —0.4 (INVl). It 
is possible to drive the equation of state closer to —1 today 
with lower values of a, although the value of A is then pushed 
to an undesirable energy range when compared with the typ- 
ical scales of particle physics. The second model INV2, which 
has Wo = —0.79 with a = 1, has been added to illustrate 
a power law potential with a dark energy equation of state 
which agrees with constraints found on wq using C MB, SN 
and large scale structure d ata (ISanchez et ^l2009t) . We also 
use the SUGRA model of iBrax fc Martini l|l999r which ex- 
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Figure 1. The dark energy equation of state as a function of 
expansion factor, 'w{a), for six quintessence models motivated by 
particle physics, which are either tracking or scaling solutions. 
The parametrization for w(a) is given in Eq. [6]and the four pa- 
rameter values which specify each model are given in Table [l] 
Note the left hand side of the x-axis is the present day. 

hibits tracking field behaviour. The potential in this case also 
contains an exponential term which pushes the dark energy 
equation of state to wo = —0.82. The 2EXP model is an ex- 
ample of a scaling solution and fea tures a double expon ential 
term in the scalar field potential ("Barreiro et al."2000|). The 
AS model suggested by Albrccht & Skordis (2000) belongs 
to the class of scaling quintessence fields. As mentioned pre- 
viously, the parameters in this potential can be adjusted to 
have the fractional dark energy density, Ode, below the nu- 
cleosynthesis bound in the early universe. The CNR model 
ijCopeland ct al. 2000) has a tracking potential where the 
scalar field rolls down to its minimum and will settle down 
to Wo = —1 after a series of small oscillations. 

Each of the quintessence models we consider is one of 
a family of such models with parameter values chosen in 
order to solve the issues of fine-tuning and coincidence, as 
well as to produce a value of uiq ~ — 1 today. These require- 
ments limit the parameter space available to a particular 
quinte ssence potential . For example, this limits the range 
of the iBrax fc Martini (1 1999) SUGRA model. The SUGRA 
model we simulate has a fixed parameter value in the super- 
gravity potential but the dark energy equation of state for 
this model does not dep end strongly on this parameter (see 
Figure 4 in lBrax fc Martin 1999. 'l. 

2.2 Parametrization of w 

Given the wide range of quintessence models in the litera- 
ture it would be a great advantage, when testing these mod- 
els, to obtain one model independent equation describing 
the evolution of the dark energy equation of state without 



having to specify the potential V{lp) directly. Throughout 
this pap er we will employ the parame trization for w pro- 
posed bv lCorasaniti fc Copel and ( 2003"), which is a general- 
isation of the method used by|Bassctt ot al. (2002.) for fitting 
dark energy models with rapid late time transitions. Using 
a parametrization for the dark energy equation of state pro- 
vides us with a model independent probe of several dark 
energy properties. The dark energy equation of state, iu(a), 
is described by its value during radiation domination, Wi, 
followed by a transition to a plateau in the matter domi- 
nated era, Wm, before making the transition to the present 
day value wo . Each of these transitions can be parametrized 
by the scale factor ttr^m at which they occur and the width 
of the transition Ar^m- 

In order to reduce this parameter space we use the 
shorter version of this parametrization for w, which is rele- 
vant as our simulations begin in the matter dominated era. 
The equation for w valid after matter-radiation equality is 

w^{a) = Wo + {uJm — Wo) X srr^ ^ 1 — -(6) 

ICorasaniti fc Copelandl l|2003h showed that this four param- 
eter fit gives an excellent match to the exact equation of 
state. Table [1] gives the best fit values for the equation of 
state parameters for the different q uintessence models taken 
from lCorasaniti fc Copelandl (|2003l ). with the addition of the 
INV2 model. The parametrization for the dark energy equa- 
tion of state is plotted in Fig. [1] for the various quintessence 
models used in this paper. 

Fig. [2] shows the evolution of the dark energy density 
with expansion factor in each quintessence model. Some of 
these models display significant levels of dark energy at high 
redshifts in contrast to a ACDM cosmology. As the AS, 
CNR, 2EXP and SUGRA models have non-negligible dark 
energy at early times, all of these could be classed as 'early 
dark energy' models. As shown in Fig.[2]both the CNR and 
the 2EXP models have high levels of dark energy at high 
redshifts compared to ACDM; after an early rapid transi- 
tion, the dark energy density evolves in the same way as 
in a ACDM cosmology. Other models, Uke the AS, INVl 
and the SUGRA models, also have non-negligible amounts 
of dark energy at early times, and after a late-time tran- 
sition, the dark energy density mimics a ACDM cosmol- 
ogy at very low redshifts. In Section 4 we will investigate 
if quintessence models which feature an early or late tran- 
sition in their equation of state, and in their dark energy 
density, can be distinguished from ACDM by examining the 
growth of large scale structure. The luminosity distance and 
Hubble parameter in the quintessence models are compared 
to ACDM in Fig. [3] and Fig. 2] respectively. In these plots 
it is clear that the CNR and the 2EXP models differ from 
ACDM only at very high redshifts. 

The adoption of a 4 variable parametrization is essen- 
tial to accurately model the expansion history over the full 
range of redshifts probed by the simulations. Using a 1 or 
2 parameter equation of state whose application is limited 
to low redshift measurements restricts the analysis of the 
properties of dark energy and cannot make use of high 
redshift m e asurem ents such as the CMB. As an example, 
ICorasanitil l|2004 ) demonstrate d that a two parameter lo g 
expansion for w{z) proposed bv lCerke fc EfstathiotJ l|2002h . 
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Table 1. The equation of state of the dark energy models s imu- 
lated, expressed in the parametrization of Corasaniti & Cope landl 
l|2003l l. The evolution of ui(a) is described by four parameters, the 
value of the equation of state today, wq , and during matter dom- 
ination era, Wm, the expansion factor, am, when the field changes 
its value during matter domination and the v^idth of the tran- 
sition. Am- We have added the INV2 model to this list as an 
example of an inverse power law potential with a value of wq 
closer to -1 than in the INVl model. 
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can only take into account a quintessence model which varies 
slowly and cann ot faithfully reprodu ce the original 'w{z) 
at high redshifts. iBassett et al.l l|2004) analysed how accu- 
rately various parametrizations could reproduce the dynam- 
ics of quintessence models. They found that parametriza- 
tions based on an expansion to first order in z or log z showed 
errors of ~ 10% at z = 1. A general prescription for w{z) 
containing more parameters than a simple 1 or 2 variable 
equation can accurately d escribe both slowly and rapidly 
varying equations of state jBassett et al."2004 l). For exam- 

J ile, th e parametrization provided bv iCorasaniti fc Copelan"dl 
20031 ) can accurately mimic the exact time behaviour of 
w{z) to < 5% for z < 10^ using a 4 parameter equation of 
state and to < 9% for z < 10"' with a 6 parameter equation. 
Finally, we note that th e para metrization for w proposed by 
ICorasaniti fc Copelandl ll2003t) is similar to the fou r param- 
eter equation of state in iLinder fc Hutered (|2005l ) (Model 
4.0) where the evolution of w is described in terms of the 
e-fold variable, A'^ = Ino, where a is the scale factor. 



2.3 The expected impact of dark energy on 
structure formation 

The growth of structure is sensitive to the amount of dark 
energy, as this changes the rate of expansion of the Universe. 
As a result, a quintessence model with a varying equation 
of state could display different large scale structure from 
a ACDM model. Varying the equation of state will result 
in different amounts of dark energy at different times. It 
has been shown that models with a larger density of dark 
energy at high redshift than ACDM have more developed 
large scale stru cture at early times, when normalised to th e 
same erg today (jCrossi fc Springeill2009l : lFrancis et al.ll2008l ). 

When the dark matter perturbations are small and the 
density contrast 5{x,t) ^ 1, the expression for the power 
spectrum as a function of time, P{k,t), is separable as 



P{k,t) 



D{to)' 



P{k,to), 



(7) 



where D{to) is the linear growth factor at the present epoch. 
The normalised gro wth factor G = DIa obe ys the following 
evolution equation l|Linder fc Jenkingl2003l ). 
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Figure 2. The dark energy density, nDE(a), as a function of ex- 
pansion factor. The INVl, SUGRA, CNR, 2EXP and AS models 
have significant levels of dark energy at early times. From z ^ 9 
until today the 2EXP and CNR models display the same en- 
ergy density as ACDM. Note the x-axis scale on this plot goes to 
z > 300 on the right hand side. 
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Figure 3. The luminosity distance in different quintessence mod- 
els compared to that in a ACDM cosmology. In this case we have 
assumed the same matter density today of f2ni = 0.26 in each of 
the models. The CNR and 2EXP models predict the same Dj^ as 
in ACDM and are overplotted. 
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Figure 4. The ratio of tlie Hubble parameter for quintessence 
cosmologies to that in ACDM. 
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where 



3 w(a) 
2 1 + X(a) 



3 l-M(a) G_ 
~ ^ 21+X{a) 'a? 



X{a) = 



fin 



-3 dlna'iu{a') 



(8) 



(9) 



and w{a) is the dynamical dark energy equation of state. The 
linear growth factor for each quintessence model is plotted in 
Fig. O In Section 14.11 we present the simulation results for 
each quintessence model where the initial conditions were 
generated using a ACDM linear theory power spectrum and 
the background cosmological parameters are the best fit val- 
ues assuming a ACDM cosmology (Stage I). The difference 
between the simulations is the result of having a different 
linear growth rate for the dark matter perturbations. 

The presence of small but appreciable amounts of dark 
energy at early times also modifies the growth rate of fluc- 
tuations from that expected in a matter dominated universe 
and hence changes the shape of the linear theory P{k) from 
the ACDM prediction. The quintessence scalar fleld can con- 
tribute at most a small fraction of the total energy density 
at early redshifts. Constraints on this amount come from big 
bang nuc leosynthesis as well as from CMB measurements. 
iBean et a l. ( 2001) found a limit of J^de < 0.045 at a ~ lO"** 
using the observed abundances of primordial nuclides and 
a constraint of f^DE < 0.39 during the rad iation domina- 
tion e ra, a ~ 10~^, from CMB anisotropies. ICaldwell et al.l 
l|2003t) discuss the parameter degeneracies which allow for 
different amounts of dark energy at early times leaving 
the position of the CMB peaks unchan ged (see Section 
14.311. Using the WMAP flrst year data, ICorasaniti et al] 
(1200J) found a limit of fioE < 0.2 at z ~ 10. Some re- 
cent parametrization dependent constraints on early dark 
energy models found the dark energy density p arameter to 
be O de < 0.02 at the last scattering surface (|Xia fc 
l2009l ). Note that all of the models we consider are consis- 
tent with this constraint, except for the AS model (see Fig. 
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Figure 5. The growth factor as a function of expansion factor. 
The upper panel shows the evolution of the linear growth factor 
in each quintessence model. In the lower panel the ratio of the 
growth factor in the quintessence models compared to ACDM is 
plotted. The growth factor in each case has been normalised to 
unity today. 



If the dark energy is not a cosmological constant, then 
there will be dark energy perturbations present, 5^ whose 
evolution will affect the dark matter power spectrum and al- 
ter the evolution equ ation in Eq. |5] (|Ferreira fc Jovcelll998l : 
IWeller fc Lewis! |2003| ). As most of the quintessence models 
we will consider display a non-negligible contribution to the 
overall density from dark energy at early times, the matter 
power spectrum is affected in two ways ([Ferreira fc Jovcd 
1 19981 : ICaldweU et aLll2003l : iDoran et al.ll2007h . In the matter 
dominated era, the growing mode solution for dark matter 
density perturbations is proportional to the expansion fac- 
tor, 5m cx a, in a universe without a scalar field component. 
In a dark energy model which has appreciable amounts of 
dark energy at early times, the dark matter growing mode 
solution on subhorizon scales is modified to become 



3m oc a 
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(10) 
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Figure 6. Linear theory power spectra at z = for dynami- 
cal dark energy quintessence models and ACDM. In this plot, 
the spectra are normalised to CMB fluctuations (on smaller 
wavenumbers than are included in the plot). The presence of a 
non-negligible dark energy density fraction at early times causes 
a scale independent suppression of growth for scales k > fccq 
where fe^q is the wavenumber corresponding to the horizon scale 
at matter radiation equality and a scale dependent suppression 
at fc < fccq. Models with high Qde at the last scattering surface 
have a lower erg today compared to ACDM if normalised to CMB 
fluctuations. 

The growth of modes on scales k > fccq, where fccq is the 
wavenumber corresponding to the horizon scale at matter 
radiation equality, is therefore suppressed relative to the 
growth expected in a ACDM universe. For fluctuations with 
wavenumbers fc < fccq during the matter dominated epoch, 
the suppression takes place after the mode enters the hori- 
zon and the growing mode is reduced relative to a model 
with fluE — 0. These t wo effects are illustrated for a scaling 
quintessence model in iFerreira fc Jovcj (|l998l l , whose Fig- 
ure 7 shows the evolution of for two wavenumbers, one 
that enters the horizon around Ocq (fc — O.lMpc"^) and one 
that comes in during the radiation era (fc = lMpc~^), in a 
universe with Que = 0.1 durning the matter dominated era. 
There is a clear suppression of growth after horizon cross- 
ing compared to a universe with no scalar field. The over- 
all result is a scale independent suppression for subhorizon 
modes, a scale dependent red tilt {us < 1) for superhorizon 
modes and an overall broading of the turnover in the power 
spectrum. This change in the shape of the turnover in the 
matter power spectrum can be clearly seen in Fig. [6] for the 
AS model. This damping of the growth after horizon crossing 
will result in a smaller erg value for the quintessence models 
com pared to ACDM i f normalised to CMB fluctuations (see 
also kunz et all l2004l l . 

We have used the publicly availabl e PPF (P arametrized 
Post-Friedmann) module for CAMB, (jPang et al. 2008). to 
generate the linear theory power spectrum. This module 
supports a time dependent dark energy equation of state 
by implementing a PPF prescription for the dark energy 
perturbations with a constant sound speed = 1. Fig. [6] 



shows the dark matter power spectra at 2; = generated 
by CAMB for each quintessence model and ACDM with 
the same cosmological parameters, an initial scalar ampli- 
tude of As — 2.14 x 10~^ and a spectral index Us = 0.96 
jSanchez et al.ll2009h . As can be seen in this plot, models 
with higher fractional energy densities at early times have a 
lower as today and a broader turnover in P{k). In Section 
14. 21 a consistent linear theory power spectrum was used for 
each quintessence model to generate the initial conditions 
for the simulations (Stage II). 

Finally, quintessence dark energy models will not nec- 
essarily agree with observational data when adopting the 
cosmological parameters derived assuming a ACDM cosmol- 
ogy. We consider how the different quintessence models af- 
fect various distance scales. We find the best fit cosmological 
parameters for each quintessence model using the observa- 
tional constraints on distances such as the measurements of 
the angular diameter distance and sound horizon at the last 
scattering surface from the cosmic microwave background. 
The method and data sets used are given in Appendix lAl and 
the corresponding simulation results which use a consistent 
linear theory power spectrum for each model together with 
the best fit cosmological parameters are presented in Section 
113] (Stage III). 



3 SIMULATION DETAILS 

We will determine the impact of quintessence dark energy 
on the growth of cosmological structures through a series 
of large N-body simulations. These simulations were car- 
ried out at the Institute of Computational Cosmology using 
a memory efficient version of the TreePM code Gadget-2, 
called L-Gadget-2 (Springel 2005). As our starting point, 
we consider a ACDM model with the following cosmological 
parameters: = 0.26, ^de = 0. 74, Ob = 0-044, h = 0.715 
and a spectral tilt of ns = 0.96 (jSanchez et al.ll2009l ). The 
linear theory rms fluctuation in spheres of radius 8 Mpc 
is set to be erg — 0.8. For each of the quintessence models the 
parametrization for the dark energy equation of state given 
in Eq. [6] was used. In the first stage we fix the cosmologi- 
cal parameters for all of the quintessence models to those 
of ACDM. As a result, some of the scalar field models do 
not match observational constraints on the sound horizon 
at last scattering or the angular diameter distance. We shall 
discuss this further in Section [4.31 using the results given in 
Appendix 1X1 

The simulations use = 646^ ~ 269 x 10*^ particles 
to represent the dark matter in a computational box of co- 
moving length 1500/i~^Mpc. We chose a comoving softening 
length of e = 50/i~^kpc. The particle mass in the simula- 
tion is 9.02 X 10^^ Mq with a mean interparticle separa- 
tion of r ~ 2.3 h~^Mpc. The initial conditions of the parti- 
cle load were set up with a glass configuration of particles. 
This arrangement is obtained by evolving a random distri- 
bution of particles with the sign of the gravitational force 
reversed (White . 1994 : Baugh et al.lll995[ l. The particles are 
perturbed from the glass using the Zeldovich approxima- 
tion which can induce small scale transients in the measured 
power spectrum. These transients die away after ~10 exp an- 
sion factors from the starting redshift l|Smith et al ] |2003l ). In 
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order to limit the effects of the initial displacement scheme 
we chose a starting redshift of 2: = 200. 

The linear theory power spectrum used to generate the 
initial conditions was created using the CAMB package of 
iLewis fc Bridia (1200^ 1. In the first stage of our calculations, 
presented in Section [4.11 the linear theory power spectrum 
used to set up the initial conditions in the quintessence mod- 
els was the same as ACDM. For the purpose of computing 
the shape of P{k) in Stage I, we have assumed that the ratio 
of dark energy density to the critical density at the last scat- 
tering surface (^jgg ~ 1000) is negligible and have ignored 
any clustering of the scalar field dark energy. In Section [4.21 
the linear theory P{k) is generated for each quintessence 
model using a modified version of CAMB which incorpo- 
rates the influence of dark energy on dark matter clustering 
at early times. In each model the power spectra at redshift 
zero have been normalised to have ag, = 0.8. Using the linear 
growth factor for each dark energy model, the linear theory 
P{k) was then evolved backwards to the starting redshift 
of z = 200 in order to generate the initial conditions for 
L-Gadget-2. Snapshot outputs of the dark matter distribu- 
tion as well as the group catalogues were made at redshifts 
5, 3, 2.5, 2, 1.5, 1, 0.75, 0.5, 0.25 and 0. The simulation code 
L-Gadget-2 has an inbuilt friends-of-friends (FOF) group 
finder which was applied to produce group catalogues of 
dark matter particles with 10 or more particles. A linking 
length of 0.2 times the mean interparticle separation was 
used in the group finder. 

We investigate gravitational collapse in the six 
quintessence models listed in Table 1 by comparing the 
evolution of the power spectrum at various redshifts. The 
power spectrum was computed by assigning the particles to 
a mesh using the cloud in cell (CIC) assignment scheme 
(|Hocknev &: EastwoodI Il98lh and then performing a fast 
Fourier transform on the density field. To restore the res- 
olution of the true density field this assignment scheme is 
corrected for by perform ing an approximate de-convolution 
jBaumgart fc Frvlll99ll ). 



4 RESULTS 

In the following sections we present the power spectrum pre- 
dictions from the three stages of simulations carried out as 
described in Section 2.3. The bottom line results are pre- 
sented in Section 4.3, in which we compare power spectra 
in ACDM with a subset of dark energy models which also 
pass the currently available observational constraints. The 
reader pressed for time may wish to skip directly to this 
section. Sections 14.11 and 14.21 show intermediate steps away 
from ACDM towards the consistent dark energy models pre- 
sented in Section [4.31 to allow us to understand the impact 
on P{k). In Section [4.11 the Friedmann equation was mod- 
ified with the quintessence model's equation of state as a 
function of redshift and a ACDM linear theory power spec- 
trum was used to generate the initial conditions for all the 
simulations (Stage I). In Section 14.21 we use a consistent 
linear theory power spectrum for each quintessence model 
(Stage II). In Section [4.31 we constrain a set of cosmological 
parameters, using CMB, BAO and SN data, for each dark 
energy model. The final stage of simulations use a consistent 
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Figure 7. Power spectra in a ACDM cosmology (orange lines) 
and AS quintessence model (green lines) at redshift , 1 an d 5. 
The red dashed lines corresponds to the [Smith et al.l (|2003[) an- 
alytical expression for the nonlinear P{k) in ACDM; blue dotted 
lines show the equivalent for the AS model. The solid black line 
is the line ar theory for ACDM at the corresponding redshift out- 
puts. The [Smith et al. I 1I2OO3II expression for the AS model has 
been scaled with the appropriate growth factor for this model at 
each redshift. 
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Figure 9. Ratio of power spectra output from the simulations 
in the six quintessence models compared to the nonlinear ACDM 
P{k) at redshift 0. Note the expanded scale on the y-axis. As 
expected, the 2EXP and CNR models show no difference from 
ACDM while the difference in the INVl, INV2, SUGRA and AS 
models is under 10% for wavenumbers k < l/iMpc~^. 
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Figure 8. The nonlinear growth of the power spectra in the various quintessence models as indicated by the key in the top left panel. 
Each panel shows a different redshift. The power spectra in each case have been divided by the ACDM power spectrum at redshift 5 
scaled to take out the difference between the ACDM growth factor at z = 5 and the redshift plotted in the panel. This removes the 
sampling variance due to the finite box size and highlights the enhanced nonlinear growth found in quintessence cosmologies compared to 
ACDM. A deviation of the power ratio from unity therefore indicates a difference in P{k) from the linear perturbation theory of ACDM. 



linear theory power spectrum for each model together with 
the best fit cosmological parameters (Stage III). 

4.1 Stage I : Changing the expansion rate of the 
Universe 

In this first stage of simulations, the same ACDM initial 
power spectrum and cosmological parameters were used for 
all models. In Fig. [7] we plot the power spectrum at red- 
shifts z = 0,1,5 in ACDM (orange lines) and in the AS 
model (green lines), together with the linear theory power 
spectra for ACDM (black lines). The AS model has a linear 
growth rate that differs from ACDM by ~ 20% at z = 5. We 
also plot the Smith ct al. (2003) 'Halofit' empirical fitting 
function for ACDM and the AS model. The Halofit function 
has been incorporated into the CAMB package and this code 
was used to generate the out put at variou s redsh ifts seen in 
Fig. [71 As this plot shows, the lSmith et al.l (l2003t ) expression 
accurately describes the evolution of the power spectrum at 
redshift in both models and at earlier times. As the nor- 
malisation and linear spectral shape is the same in these two 
models, Halofit accurately reproduces the nonlinear power 
in each model at various redshifts once the appropriate lin- 



ear growth factor for the dark energy model at that redshift 
is used. The Smith et al. expression agrees with the sim- 
ulation output at z = to within 4% for k < l/iMpc~^ 
for both the quintessence model and ACDM. At higher red- 
shifts, the difference between the simulation output and the 
Halofit prediction for all the models is just under 10% on 
scales k < 0.3/iMpc~^ at z = 5. 

To highlight the differences in the power between the 
different models, we plot in Fig. [8] the measured power 
divided by the power at z = 5, after scaling to take 
into account the difference in the linear theory growth 
factors for the output redshift and z = 5, for ACDM. 
Thi s removes the sampl i ng va riance from the plotted ra- 
tio (|Baugh fc Efstathiou|[l994h . A ratio of unity in Fig. \E\ 
would indicate linear growth at the same rate as expected 
in ACDM. 

Fig. [8] shows four epochs in the evolution of the power 
spectrum for all of the quintessence models and ACDM. The 
black line in the plot shows the P{k) ratio for ACDM (note 
the yellow curve for the CNR model is overplotted) . Non- 
linear growth can be seen as an increase in the power ratio on 
small scales, k > 0.3/iMpc~^ at z = 3 and k > O.lftMpc"^ 
at z = 0. Four of the quintessence models (INVl, INV2, 
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Figure 10. The ratio of the quintessence model power spectra to the ACDM power spectrum output from the simulations at three 
values of the linear growth factor D = I, D = 0.5 and D = 0.3. Each panel shows the results of this exercise for the AS, CNR, 2EXP 
and SUGRA quintessence models. The growth factors correspond to z = 3.4 (D = 0.3), z = 1.6 {D = 0.5) and z = {D = 1) for ACDM. 
For each model, the choice of growth factor corresponds to slightly different redshifts, with the biggest difference being for the INVl 
model. A ratio of unity would indicate that the growth factor is the only ingredient needed to predict the power spectrum in the different 
quintessence models. Note the expanded scale on the y axis. 



SUGRA and AS) differ significantly from ACDM for z > 0. 
These models show advanced structure formation i.e. more 
power than ACDM, and a large increase in the amount of 
nonlinear growth. All models are normalised to have as = 
0.8 today and as a result all the power spectra are very 
similar at redshift zero in Fig. [8] There are actually small 
differences between the quintessence models at z = as seen 
on the expanded scale in Fig. [51 This increase in nonlinear 
power at small scales in the quintessence models is due to 
the different growth histories. 

The power spectra predicted in the 2EXP and CNR 
models show minor departures from that in the ACDM cos- 
mology. This is expected as Figs. [T] and [2] show the equations 
of state and the dark energy densities in these two models 
are the same as ACDM at low redshifts and all three sim- 
ulations began from identical initial conditions. It could be 
possible to distinguish these two models from the concor- 
dance cosmology at higher redshifts if we do not ignore the 
dark energy perturbations or changes in the growth factor 
which alter the form of the linear theory power spectrum. 
We shall discuss this more in the next stage of our simula- 
tions in Section [4.21 

Finally, we investigate if the enhanced growth in the 



power spectrum seen in Fig. [8] in the quintessence models 
is due solely to the different linear growth rates at a given 
redshift in the models. In order to test this idea, the power 
spectrum in a quintessence model and ACDM are compared 
not at the same redshift but at the same linear growth factor 
0. As the growth rates in some of the quintessence models are 
very different from that in the standard ACDM cosmology, 
the power spectra required from the simulation will be at 
different output redshift in this comparison. For example, 
the normalised linear growth factor is D = 0.5 at a redshift 
of z — 1.58 in a ACDM model and has the same value at 
2 = 1.82 in the SUGRA model, at z = 1.75 in the AS 
model and at z — 2.25 in the INVl quintessence model. In 
Fig. [TO] we show the power spectrum of simulation outputs 
from the INVl, AS, SUGRA and CNR models divided by 
the power spectrum output in ACDM at the same linear 
growth rate. We ran the simulations taking three additional 
redshift outputs where the linear growth rate had values 
of D = 1,D = 0.5 and D = 0.3. It is clear from Fig. [10] 
that scaling the power spectrum in this way can explain the 
enhanced linear and most of the excess nonlinear growth 

^ We thank S. D. M. White for this suggestion. 
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seen in Fig.[8]for scales k < 0.1/iMpc ^ . For example, in the 
INVl model the enhanced nonlinear growth, on scales k ~ 
0.3/iMpc"^ at fixed D = 0.3, difTers from ACDM by at most 
5% in Fig. [To] as opposed to at most 30% at z = 5 in Fig. [S] 
At earlier redshifts when the linear growth rate is D — 0.3, 
the nonlinear growth in the quintessence models agrees with 
ACDM on smaller wavenumbers k < 0.3/iMpc~^. As in Fig. 
[S] the CNR model shows no difference from ACDM when 
plotted in this way. 

Note in Fig. [10] the INVl model has less nonlinear 
growth at _D = 0.3 and D = 0.5 compared to the AS model. 
The AS and SUGRA models have a growth rate of D=0.5 
at lower redshifts compared to the INVl model and so are 
at a later stage in their growth history. The INVl model has 
a growth rate of D = 0.5 a.t z — 2.25 whereas for the AS 
model this occurs at z = 1.75 and at 1.82 for the SUGRA 
model. The reason for the success of this simple model - 
matching the growth factor to predict the clustering - can 
be traced to the universality of the mass function, which we 
discuss in Section [4.41 In this Stage I calculation, the mod- 
els have the same mass function when plotted at the epoch 
corresponding to a common growth factor. This means that 
the two-halo contribution to the clustering is therefore the 
same. Can this simple halo picture of the clustering also 
explain the clustering on small scales (high A;)? Although 
the abundance of haloes in the models is the same at the 
epochs corresponding to a given value of the growth factor, 
the concentrations of the haloes will not be the same. In cos- 
mologies where the haloes formed at a higher redshift (i.e. 
roughly the redshift corresponding to a particular value of 
D), one would expect these haloes to have higher concentra- 
tions than their counterparts in the other models ijEke et al.l 
I2OOII ). A higher concentration would be expected to yield 
stronger nonlinear clustering and hence more power at high 
k in Fig. [10] Unfortunately our simulations do not have the 
resolution to probe the required range of wavenumbers to 
uncover this behaviour. The ratios plotted in Fig. [10] stop 
at wavenumbers approximately equivalent to the collapsed 
radius of a massive halo. 

Hence, it seems that scaling the power spectrum us- 
ing the linear growth rate can be used to predict the linear 
growth in the quintessence dark energy simulations and can 
reproduce some of the nonlinear growth at early redshifts. 
In Fig. [TO] there are still some differences in the small scale 
growth in quintessence models compared to ACDM which 
cannot be explained by the different linear growth rates. We 
find that nonlinear evolution is not just a function of the 
current value of the linear growth rate but also depends on 
its history through the evolution of the coupling between 
long and short-wavelength modes. 



4.2 Stage II: Use of a self-consistent linear theory 

P{k) 

We have run the simulations presented in the previous sec- 
tion again but this time using the appropriate linear theory 
P{k) for each model (shown in Fig.[S]) normalised to as = 0.8 
today (Stage II). After normalising the power spectra in this 
way, the difference between the quintessence models P{k) 
and ACDM can be seen in Fig. [TT] The INV2 model was 
not included in this set of simulations as there is a negli- 
gible difference in the linear theory power spectrum from 
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Figure 11. Ratio of linear theory power spectra for quintessence 
models shown in Fig. [6]to that in ACDM. In this plot each P{k) 
has been normalised so that ag = 0.8 today; this is the normali- 
sation used in our simulations. 



ACDM. Note [Francis etld] (|2008l ') also generate the linear 
theory power spectrum for 'early dark energy ' models and 
norma lise all P{k) to have the same erg today. i Francis et al.l 
l|2008l ) make an equivalent plot to Fig. [TT] but find a de- 
crease in this ratio with decreasing scale (fc > 0.2/iMpc~^), 
using the paramet r izatio n for early dark energy proposed by 
[Poran fc Robber^ (|2006h . in contrast to the ratio of unity we 
find on small scales in Fig. 1111 This difference is due to the 
different parametrizations used for the dark energy equation 
of state, as a ratio of unity is obtained on small scales for the 
same 'early dark energy' mo del using the parametrization 
sugg ested bv lWetterichfl|2004l ') (M. Francis, private commu- 
nication) . 

In the first row of Fig.[T2]we plot the power spectrum for 
the Stage II SUGRA model at z — 0, 1, and 3 divided by the 
simulation output in ACDM at 2: = 5 as in Fig. [8] (red dashed 
lines). The result from Fig. [8] Stage I SUGRA, is also plot- 
ted here to highlight how changing the spectral shape affects 
the nonlinear growth in the simulations. On large scales the 
growth is not modified by the altered spectral shape. The 
growth of perturbations on small scales in the simulation is 
affected by the modified linear theory used in the initial con- 
ditions. Normalising the power spectra to erg = 0.8 results in 
more power on large scales in the quintessence models com- 
pared to ACDM, as can be seen in Fig. 1111 This enhanced 
large scale power couples to the power on smaller scales and 
results in a small increase in the nonlinear power spectrum 
for k > 0.1/iMpc~^ in the Stage II SUGRA simulation com- 
pared to the one using ACDM linear theory P(k) in Stage 
I. 

In the second row of Fig. [12] we plot the power spectrum 
for the Stage II AS model as green dashed lines at 2 = 0, 1, 
and 3 divided by the simulation output in ACDM at 2 = 5 as 
in Fig. [8] The growth of dark matter perturbations is greatly 
suppressed in the AS model due to the large fractional dark 
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Figure 12. Ratios of power spectra for the SUGRA (first row), AS (second row) and CNR (tliird row) quintessence model compared to 
ACDM from tlie 3 stages of simulations in this paper. The plot shows the growth in the quintessence models using ACDM linear theory 
P{k) in the initial conditions in black (Stage I) and using a self consistent linear theory P{k) for each quintessence model (dashed colored 
line) (Stage II). The dotted lines shows the P(k) ratio from the simulation for the quintessence models using the best fit parameters in 
Table A3 (Stage III). The power spectra in each case have been divided by the ACDM power spectrum at redshift 5 with appropriate 
scaling of ACDM growth factors. The linear theory power spectra in each case has been normalised to erg = 0.8. 



energy density at high redshifts. After fixing erg = 0.8, there 
is more power on large scales in the AS model compared to 
ACDM. As in the first row of Fig. [T^lthere is a small increase 
in nonlinear power for the AS model in Stage II. Although 
the excess large scale power is significantly larger than in 
the SUGRA model case, it does not result in more nonlinear 
power on small scales through mode coupling, as can be seen 
in the panels in the second row in Fig. 1121 The linear theory 
power spectrum for these quintessence models has a scale 
dependent red tilt on large scales which shifts the position of 
the BAO peaks which is the origin of the oscillation apparent 
in the second row of Fig.[T2]at z = 3. The difference in BAO 
peak positions is very prominent when we plot the ratio of 



the power spectrum in the AS model to the ACDM power 
spectrum and can be clearly seen in Fig. [121 



4.3 Stage III: Consistency with observational 
data 

In this section we present the power spectra results in 
ACDM and a subset of the dark energy models, measured 
from simulations which use a consistent linear theory power 
spectrum for each model together with the best fit cosmo- 
logical parameters. We have simulated the SUGRA, AS and 
CNR models using the best fit cosmological parameters from 
Table A3 and the linear theory power spectrum specific to 
each model as discussed in Section 12.31 We chose to simu- 
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late these three models following the analysis and results of 
Sections 14.11 14.21 and Appendix |X] Any of the dark energy 
models listed in Section 2.2 which showed similar results in 
Section [4.21 to ACDM and similar cosmological parameters 
in Appendix 1X1 have not been simulated again. 

Table A3 in Appendix shows the best fit values for 
Qmh^, f2b/i^ and Ho for each quintessence model, found by 
minimising xLai = Xwmap+sn+bao- The SUGRA, AS and 
CNR models had the biggest improvement in the agreement 
with observational constraints, on allowing Qmh^ , 0,hh^ and 
Ho to vary. The results for the SUGRA, AS and the CNR 
model are shown as dotted coloured lines in Fig. [12] and are 
referred to as Stage III in the legend to distinguish them 
from the results of Sections 14.11 and 14.21 which are also plot- 
ted. In each row we show the simulation outputs at 2; = 0, 1 
and 3. The simulation results for each quintessence model 
uses the models linear theory and the best fit parameters 
from Table A3. Using the best fit parameters for each model 
together with the correct linear theory changes the growth 
of structure in the simulation. 

In Fig.[T2]the measured power spectrum for each model 
is divided by the power for ACDM at 2; = 5 which has been 
scaled using the difi'erence in the linear growth factor be- 
tween z = 5 and the redshift shown. Plotting the ratio in 
this way highlights the differences in growth between the 
quintessence models and ACDM as well as removing sam- 
pling variance. 

The measured power for the SUGRA model is plotted 
in the first row in Fig. 1121 The power spectra have all been 
normalised to erg = 0.8 resulting in a large increase in the 
large scale power {k < 0.1/iMpc~^) seen in Fig. [12] com- 
pared to ACDM. There is a large increase in the linear and 
nonlinear growth in this model at 2 > (dotted red line) 
compared to ACDM (dot-dashed grey line). The second row 
in Fig [T^ shows there is a significant enhancement in the 
growth in the AS power spectrum measured compared to 
ACDM for 2 < 3. The power measured from the simulations 
of the CNR model are plotted in the third row of Fig. 1121 
We find there is a small reduction in the amount of linear 
and nonlinear growth in this model compared to ACDM. 

In Fig. [T2] we also plot the simulation results for these 
three models from Section 4.1 (Stage I), where ACDM lin- 
ear theory was used in the initial conditions, (black lines). 
The dashed coloured lines show the simulation results from 
Section 4.2 (Stage II), where the quintessence model linear 
theory was used. The SUGRA power spectrum measured 
in Stage III has less nonlinear growth at high redshifts com- 
pared to the SUGRA P{k) from Stage I or II due to changes 
in the spectral shape. The measured power for the AS model 
using the best fit parameters (Stage III) shows enhanced 
growth on all scales compared to the power for the AS model 
in Stage I (using ACDM parameters and linear theory P{k)) 
or Stage II (using ACDM parameters). 

These results show the importance of each of the three 
stages in building up a complete picture of a quintessence 
dark energy model. Models whose equation of state is very 
different from ACDM at low redshifts, for example the 
SUGRA and the AS model, show enhanced nonlinear growth 
today compared to ACDM. Models whose equation of state 
is very different to ACDM only at early times, for exam- 
ple the CNR model, will show no difference in the nonlin- 
ear growth of structure if we use the ACDM spectral shape 



(Stage I). In Stage II and III the shape of the power spec- 
trum in the CNR model has changed and is very different 
to ACDM on large scales as can be seen in Fig. 1121 Using 
the best fit cosmological parameters for this model we find 
a very small reduction (< 2%) in the nonlinear growth at 
2 = compared to ACDM. 

4.4 Mass function of dark matter haloes 

In this section we present the mass function of dark matter 
haloes in the quintessence models using the three stages of 
sim ulations discussed i n Secti ons 14.11 14.21 and 14.31 

jPress fc Schechte"rl l| 19741 ) (hereafter P-S) proposed an 
analytical expression for the abundance of collapsed objects 
with mass M in the range M to M + dM at redshift 2, based 
on the spherical collapse model in which a perturbation can 
be associated with a virialised object at 2 = 2', if its density 
contrast, extrapolated to 2 = 2' using linear theory, exceeds 
some threshold value, Sc, the critical linear density contrast. 
It has been shown that the P-S approach fails to reproduce 
the abundance of haloes found in simulations, overpredict- 
ing the number of haloes below the characteristic mass M* 
and underpredicting the abundance in the high mass tail 
'^■^ ' • " Rees '1988'; 'W hite et al.1 Il993l : iLacev fc Coll 



(Efetathii . 

fl994; Eke ct al. 1996; Govcr nato et al 



|199S ) 



It is thought that the main cause of this discrepancy is 
the spherical collapse approximation, as the perturbations 
in the density field are inherently triaxial. After turnaround, 
each axis may evolve separately until the final axi s col- 
lapses and the o bject virialises. ISheth et al.l (|200ll ) and 
ISheth fc TormenI ||2002D (hereaft^er S-T) modified the P-S 
formalism, replacing the spherical collapse model with ellip- 
soidal collapse, in which the surrounding shear field as well 
as the i nitial overdensity d etermines the collapse time of an 
object. ISheth et al. I (I2OOII) found a universal mass function 
for any CDM modeh ljenkins et al.l (|200lf ) found a univer- 
sal empirical fit to the form of the mass function measured 
from a suite of cosmological simulations. The Jenkins et al. 
mass function can accurately predict halo ab undances over 
a range of cosmologies and redshifts (see also 1 Warren et al.l 
l2006l : lReed et al.ll2007l and lCrocce et al.ll2009l ). 

We use a friends-of-friends (FOF) halo finder, with a 
constant linking length of b = 0.2, to identify haloes in all 
cosmologies. In Fig. [13] we plot groups containing 20 par- 
ticles or more to ensure that the systematic uncertainties 
in the mass function are at or below the 10% level; tests 
show t hat 90% or more of such haloes are gravitationally 
bound ISpringel et al.ll2005D . The first row in Fig. 1 131 shows 
the mass function for SUGRA and ACDM at 2 = 0, 1 and 
2. The filled red squares represent the mass function from 
Stage III of the simulations where a consistent linear the- 
ory and cosmological parameters were used for the SUGRA 
model. The mass function for ACDM (open black circles) 
and the SUGRA model are plotted together with the Jenk- 
ins et al. mass function shown in black (red) for ACDM 
(SUGRA). The S-T mass function is shown in the top left 
panel in the first row of this figure (blue dashed line) for 
comparison. The abundances in both ACDM and SUGRA 
agree with each other at redshift and with the Jenkins 
et al. and S-T models, although the fitting formulae seem 
to slightly under-predict the number of haloes at the high 
mass end (M > lO^^/i"^ Mq). In the first row of Fig. [l3 
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Figure 13. Dark matter halo mass functions for the SUGRA (first row) and AS (second row) quintessence models compared with that 
in ACDM from the Stage III simulations at 2 = 0, 1 and 2. The mass function in ACDM is shown as open black circles throughout 
this plot. In the first row the red filled squares show the mass function from the simulation for the SUGRA model using the best fit 
parameters in Table A3 (Stage III). Underneath each panel in the first row we plot the log of the ratio between the measured mass 
function for ACDM (open black circles) and Stage III SUGRA (red squares) and the Jenkins mass function for ACDM. In the second row 
the green filled squares show the mass function from the simulation for the AS model using the best fit parameters in Table A3 (Stage 
III). For the AS Stage III simulation, VL^h? = 0.086, giving rise to a change in the spectral shape of the linear theory power spectrum. 
As a result, there are fewer low mass halos and a similar number of high mass haloes at 2 = compared to ACDM {Q,uih? = 0.1334). 
The difference between the Jenkins et al. mass function for ACDM and the measured mass function for ACDM (open black circles) and 
Stage III AS (green squares) is plotted underneath each panel in the second row. The black horizontal line indicates a ratio of unity in 
the ratio plots. In the first and second rows the solid black (red/green) lines are the predicted abundances in the ACDM (SUGRA/AS) 
model using the Jenkins et al. fitting function at various redshifts. In the top left panel, for reference, we have also plotted the Sheth & 
Tormen mass function (blue dashed line) for ACDM. 



the number of haloes in the two models start to differ at 
z = 1, and at z = 2 there is a large difference in the mass 
functions. The linear growth factor for the SUGRA model 
together with the best fit cosmological parameters from Ta- 
ble A3 have been used to obtain the Jenkins et al. fit at the 
earlier redshifts. The Jenkins et al. fit describes the data 
slightly better at the high mass end at higher redshifts than 
the S-T prescription. This is as expected as the Jenkins et 
al. fit was explicitly tested at the high mass end of the mass 
function. Each model shows only small (< 20%) differences 
between the measured value and the Jenkins et al. fitting 
formula for M < IQ^^H'^Mq at 2 = 0. Underneath each 
panel in the first row in Fig. 1131 we plot the ratio between 
the measured mass function for ACDM and the SUGRA 
model in Stage III, and the Jenkins at al. mass function for 
ACDM. 



The second row of Fig. [13] repeats this comparison for 
the AS model. In this row the mass function for ACDM 
(open black circles) and the AS model from Stage III (green 
squares) of the simulations &t z = 0, 1 and 2 are plotted. 
The Jenkins et al. mass function for ACDM (black line) and 
the AS model for Stage III (green line) are also plotted. The 
AS model has a greater abundance of halos than ACDM 
at z = 2. For the Stage III simulation, the AS model has 

fim/l^ ^ Q Qgg 

giving rise to a change in the spectral shape 
of the linear theory power spectrum from ACDM linear the- 
ory (ilmft^ ~ 0.133). As a result there are fewer low mass 
halos and a similar number of high mass haloes at z = 
compared to ACDM. This change accounts for the decrease 
in the mass function for M < IQ^^h'^MQ seen at 2 = 
in the AS model (green squares). At z = 0, there are only 
small (< 20%) differences between the measured value and 
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Figure 14. The halo mass function for the SUGRA and AS model and ACDM at 2 = and 1 compared to the I Jenkins et al] ll200lh 
analytic fit. The Jenkins et al. mass function is plotted as solid black (red/green) lines for ACDM (SUGRA/AS). Underneath each 
panel the ratio of the mass function measured from the simulation and the Jenkins et al. mass function is plotted for all models. Note a 
logarithmic scale is used on the y axis in the ratio plots. 



the Jenkins et al. fitting formula for M < IO^^H'^Mq for 
ACDM and the AS model from Stage III. The ratio between 
the Jenkins et al. mass function for ACDM and the measured 
mass function for ACDM and the AS model from Stage III 
is plotted underneath each panel in the second row in Fig. 

M 

Only the SUGRA and AS models are plotted in Fig. [13] 
but similar differences in halo abundances are seen in the 
INV models compared to ACDM, whilst only negligible dif- 
ferences with A CDM were found in the mass functions of 
2EXP and CNR. ICrossi fc Springell ()2009( ) found similar re- 
sults for the mass function over the range 
in an 'early dark energy' model, using much smaller volume 
simulations than ours. They found a higher number den- 
sity of haloes corresponding to groups and clusters in non- 
standard dark energy models at high redshifts compared to 
ACDM, while at z = the models all agreed with one an- 
other. We find similar results although using the cosmolog- 
ical parameters from Table A3 for each quintessence model 
can give different abundances at z = in those models com- 
pared to ACDM because although as is the same the shape 
of the linear theory can be different. Also, we have been able 
to probe a higher mass range for the dark matter haloes. 
The high mass end of the mass function is very sensitive to 
changes in the current value of the linear growth factor in 
the different cosmologies. 

In Fig.[T3]we plot the fraction of the total mass in haloes 
of mass M rather than simply the abundance as shown in 
Fig. 1131 We compare the Jenkins et al. analytic fit to our 
simulated halo mass functions in the SUGRA and AS models 
and in ACDM at 2: = and 1 in Fig. [HI In this plot the 
quantity lna~^(M, 2) is used as the mass variable instead of 
M, where a^{M, z) is the variance of the linear density field 
at 2 = 0. This variance can be expressed as 

a^{M,z) = ^^ r eP{k)W\k-M)Ak, (11) 
Jo 

where W{k; M) is a top hat window function enclosing a 



mass M, D{z) is the linear growth factor of perturbations 
at redshift 2 and P{k) is the power spectrum of the linear 
density field. Plotting different masses at different redshifts 
in this way takes out the redshift dependence in the power 
spectrum. Note a large value of \na~^{M, z) corresponds to 
a rare halo. Using this variable, Jenkins et al. found that 
the mass function at different epochs has a universal form, 
for a fixed power spectrum shape. Note that in our case, the 
Stage III simulations have somewhat different power spec- 
tra, which account for the bulk of the dispersion between the 
simulation results at the rare object end of Fig. [TJ] in Stage 
I, the simulation results agree with the Jenkins et al. univer- 
sal form to within 25% at lncr~^ = 1.0. As shown in Fig. 1141 
we find the Jenkins et al. fitting formula is accurate to ~ 20% 
at 2 = for all the models in the range M < 10^^/i"^Mq. At 
higher redshifts the measured mass function for the SUGRA 
model and ACDM differ from the Jenkins et al. mass func- 
tion by ~ 30% over the same mass range while for the AS 
model the difference is ~ 50% at 2 = 1. In previous work, 
iLinder fc Jenkitii (|2003l ) also found that the predicted mass 
function for a SUGRA-QCDM simulation, which would be 
the equivalent of our Stage I simulations, was well fit (within 
20%) by the Jenkins et al. formula. 

4.5 The appearance of baryonic acoustic 
oscillations in quintessence models 

In this section we examine the baryonic acoustic oscillation 
signal in the matter power spectrum for the AS, SUGRA 
and CNR models. IXngulo et al.l (|2008l ) presented a detailed 
set of predictions for the appearance of the BAO signal in 
the ACDM model, covering the impact of nonlinear growth, 
peculiar velocities and scale dependent redshift space dis- 
tortions and galaxy bias. Here we focus on the first of these 
effects and show power spectra in real space for the dark 
matter. We do not consider the INVl model as it is not con- 
sistent with observational constraints (Appendix A), or the 
INV2 or 2EXP models as they are indistinguishable from 
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Figure 15. The ratio of tiie distance measure Dv{z) (left panel) and tlie ratio of rs{zd)/ Dv (right panel) for fou r quintessence models 
compared to ACDM as indicated by the key in the right hand panel. The grey circles are estimate points from iPercival et alj (|2007| ) 
at z = 0.2 and z = 0.35 measured using the observed scale of BAO calculated from the SDSS and 2dFGRS main galaxy samples. 
ISanchez et al. I ll2009h combined CMB data with information on the shape of the redshift space correlation function using a larger LRG 
dataset and fo und D^{z = 0.35) = 1300 ± 31 Mpc and rg(z^) / Dy = 0.1185 it 0.0032 at z = 0.35 (blue squares). The data points from 
iPercival et al] (i2009.) for D„ and rs{zd)/Dy at z = 0.275 using WMAP 5 year data + SDSS DR7 are plotted as black triangles. 



ACDM, and hence were not simulated again in Stage III 
(Section 1331). 

In Stage I of our simulations (Section l4.ip . we would ex- 
pect the linear theory comoving BAO for the quintessence 
models to be identical to ACDM as the same linear theory 
power was used for all models. In Stage II (Section 14. 2p . 
some of the quintessence models have large amounts of dark 
energy at early times which will alter the sound horizon in 
these models compared to ACDM (see Table A3), and as 
a result we would expect to see a corresponding shift in 
the BAO peak positions. The best fit cosmological parame- 
ters found in Stage III were derived using CMB, BAO and 
SN distance measurements (see Appendix Stage III of 
our simulations (Section 14. 3p uses these parameters and we 
would expect models with the same BAO distance measures 
to have the same peak pattern in the matter power spectrum 
as ACDM. 

The baryonic acoustic oscillations are approximately a 
standard r uler and depend on t he sound horizon, r^, given 
in Eq. IA3I (jSanchez et al.l |2008| ) . The apparent size of the 
BAO scale depends on the distance to the redshift of obser- 
vation and on the ratio Vs/ D^, where _D„ is an effective dis- 
tance measure which is a combination of Da and H, given in 
Eq. lA6l In most quintessence models, remains unchanged 
unless there is appreciable dark energy at last scattering. 
Models which have the same ratio of Vs/Dv are impossible 
to distinguish using BAO. 

To calculate the power spectrum for a galaxy red- 
shift survey, the measured angular and radial separations of 
galaxies pairs are converted to co-moving separations and 
scales. This conversion is dependent on the cosmological 



model assumed in the analysis. These changes can be com- 
bined into the single effective measure, _D„. Once the power 
spectrum is calculated in one model we can simply re-scale 
P{k) using Dv to obtain the power spectrum and BAO peak 
posit ions in another cosmological model (see lSanchez et aU 
boogh . In the left panel of Fig. 1151 we plot the ratio of D^ in 
four q uintessence models compared to ACDM up to z = 1.5. 
IPerci val ct al. (2007. ) found Z)„ = 564±23/i"^Mpc at z = 0.2 
and Dy = 1019 ± 42/i"^Mpc at z = 0.35 using the observed 
scale of BAO measured from the SDSS DR5 galaxy sample 
and 2dFGRS. These data points are plotted as grey circles 
in Fig. 1151 Note that at face v alue none of th e models we 
consider are consistent with the IPercival et ah I (|2007, ) point 
a,t z = 0.35. These authors report a 2.4(t discrepancy be- 
tween their results using BAO and the constraints available 
at the time from supernovae. The blue square plotted in the 
left panel in Fig. 1151 i s the constraint _D„ = 1300 ± 31 Mpc 
at z = 0.35 found bv lSanchez et all (|2009l ). This constraint 
was found using a much larger LRG dataset and improved 
modelling of the co rrelation function on large scales. The 
constraint found by I Sanchez et al.l (|2009l ) using CMB and 
BAO data is fully consisteii t with CMB and SN results. The 
results from IPercival et al.l (|2009l ) for 7?„ and rs{zd)/Dv at 
z = 0.275 using WMAP 5 year data together with the SDSS 
data releaseJ7 galaxy sample are also plotted (black trian- 
gles). The Percival et al. (20091) results are in much better 
agreement with those of lSanchez et al.l (|2009l ). 

Over the range of redshifts plotted in Fig. [15] the dis- 
tance measure, in the AS, 2EXP and CNR models difi'er 
from ACDM by at most 2% and is < 1% in these models 
for z < 0.2. Re-scaling the power spectrum for these dark 
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Figure 16. The real space power spectrum for the AS model on large scales at 2 = (left) and 2 = 3 (right). All power spectra have 
been divided by a smoothed linear 'no-wiggle' theory P{k) for ACDM. The factor, /, removes the scatter of the power measured in the 
simulation around the expected linear theory power. Stage I in our simulation is represented by grey circles, Stage II is represented by 
open blue squares and Stage III results are shown as green triangles. The black solid line represents the linear theory power spectrum in 
ACDM divided by the smooth reference spectrum. The vertical dashed (dotted) lines show the position of the first two acoustic peaks 
(positions ±5%) for a ACDM cosmology. 
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Figure 17. The real space power spectrum for the SUGRA model on large scales at 2 = and 2 = 3. All power spectra have been 
divided by a smoothed linear theory P{k) for ACDM. Stage I in our simulation is represented by grey circles. Stage II is represented by 
open purple squares and Stage III results are shown as red triangles. The black solid line represent the linear theory power spectrum in 
ACDM divided by the smooth reference spectrum. The vertical dashed (dotted) lines show the position of the first two acoustic peaks 
(positions ±5%) for a ACDM cosmology. 
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energy cosmologies would result in a small shift ~ 1% in the 
position of the peaks at low redshifts. The value of Z)„ in 
the SUGRA model differs from ACDM by at most 9% up 
to z = 1.5. The right panel in Fig. 1151 shows the ratio of 
rs{zd)/Dv in the quintessence models compared to ACDM, 
where rs is the co-moving sound horizon scale at the drag 
redshift, Zd, which we discuss in Appendix|X] The value of 
ra{zd)/Dv can be constrained using the position of the BAO 
in the power spectrum. In t he right panel of Fig. TSJthe grey 



symbols are the results from lPercival et ahl l|2007 ) at z = 0.2 
and z = 0.35. From this plot it is clear that the SUGRA and 
AS model are within the la limits at z = 0.2. The 2EXP 
and CNR model lie just outside the la errors at z = 0.35. 
Note the value of rs(zd)I Dv for ACDM a.t z — 0.35 also lie 
outside the la errors, see iPercival et al] (|2009l ) for more de- 
tail. The blue square plotted in the right panel in Fig. llSl is 
rs{zd)IDy = 0.1185 ± 0.0032 at z = 0.35 and was obtained 
using information on the redshift space correla tion function 
together with CMB data (|Sanchez et al.ll2009l '). 

In Fig. [16] and [17] we plot the 2 = and z — 3 
power spectra in the AS and SUGRA models divided by 
a linear theory ACDM reference spectrum which has been 
smoothed using the coarse rebinning method prop o sed b y 
iPercival et al! }20o3) and refined by lAngulo et all (|2008D . 
After dividing by this smoothed power spectrum, the acous- 
tic peaks are more visible in the quasi-linear regime. In 
Figs. [TB] and 1171 the measured power in each bin has been 
multiplied by a factor, /, to remove the scatter due to 
the small number of large scale modes in the simulation 
l|Baugh fc Efstathioul [iQM : ISpringel et al.ll2005l ). This fac- 
tor, / = P(A;)ii„car/f'(fc)N-body, Is thc ratlo of the expected 
linear theory power and the measured power in each bin 
at 2 = 5, at which time the power on these scales is still 
expected to be linear. Multiplying by this correction fac- 
tor allows us to see the onset of nonlinear growth around 
k ~ 0.15/iMpc~^ more clearly. 

In Fig. [16] we plot the AS (SUGRA) power spec- 
trum as grey circles from Stage I, blue (purple) squares from 
Stage II and green (red) triangles from Stage III. The black 
line represents the linear theory power in ACDM divided 
by the smooth reference spectrum. In both plots and for all 
power spectra, the same reference spectrum is used. The ref- 
erence is a simple 'wiggle-free' CDM spectru m, with a form 
controlled by the shape parameter F = f2mft (|Bardeen et al.l 
Ll986). The difference between the AS and ACDM hnear the- 
ory, as shown in Fig. 1111 results in an increase in large scale 
power on scales k < 0.04/iMpc~^. The vertical dashed (dot- 
ted) lines show the first two positions of the acoustic peaks 
(positions ±5%) for a ACDM cosmology. 

As shown in Fig. 1161 we find that the position of the first 
acoustic peak in the AS model from Stage I is the same as 
in ACDM. The position of the first peak for the AS model, 
measured in Stage II of our simulations (blue squares), is 
slightly shifted (~ 4%) to smaller scales compared to ACDM 
as the sound horizon is altered in the AS model. In Stage III, 
when the best fit cosmological parameters for the AS model 
are used, the sound horizon in the AS model and ACDM are 
very similar at 2 ~ 1090 and there is a very small (< 1%) 
shift in the position of the first peak (green triangles). As 
there is less nonlinear growth at 2 = 3 the higher order peaks 
are more visible in the right-hand plot in Fig. 1161 

In Fig. 1171 the SUGRA power spectrum from Stage I, II 



and III are plotted. The SUGRA P{k) from Stages I and II 
have identical peak positions to ACDM as the sound horizon 
is the same as in ACDM in these cases. There is a shift 
(~ 5%) in the position of the first peak in the SUGRA model 
using the P{k) measured in Stage III. Note the units on the x 
axis are h/Mpc and from Table A3, h = 0.67 for the Stage III 
SUGRA model compared to h = 0.715 for ACDM. On smaU 
scales the BAO signature is damped due to more nonlinear 
structure formation at 2 = compared to 2 = 3 as shown in 
Fig 1171 We find a large increase in the power in the region 
of the second peak, k ~ 0.15/iMpc~^ in both the AS and 
SUGRA models, measured in Stage III, compared to ACDM. 
For brevity we have not included the plots of the power 
spectra for the CNR model showing the baryonic acoustic 
oscillations. We find identical peak positions in ACDM and 
this model in all stages at 2 = 0. 

The AS and SUGRA model are very different to ACDM 
at late times and as result they affect the growth of structure 
at 2 > as seen in Sections 14.31 and 14.41 We have found that 
models like this do not necessarily have different BAO peak 
positions to ACDM in the matter power spectrum. These re- 
sults suggest that distinguishing a quintessence model, like 
the AS model used in this paper, using measurements of 
the BAO peak positions in future galaxy surverys, will be 
extremely difficult. The BAO peak positions for the CNR 
model will be shifted by at most 2% in the range 2 < 1.5 
compared to ACDM after re-scaling the power spectra by 
_D„. In conclusion it is possible to have quintessence cos- 
mologies with higher levels of dark energy at early times 
than in ACDM and still measure the same peak positions 
for the BAO in the matter power spectrum. 



5 CONCLUSIONS AND SUMMARY 

Observing the dynamics of dark energy is the central goal 
of future galaxy surveys and would distinguish a cosmolog- 
ical constant from a dynamical quintessence model. Using 
a broad range of quintessence models, with either a slowly 
or rapidy varying equation of state, we have analysed the 
infiuence of dynamical dark energy on structure formation 
using N-body simulations. 

We have considered a range of quintessence models that 
can be classified as either 'tracking' models, for example the 
SUGRA and INV models, or 'scaling' solutions, such as the 
AS, CNR or 2EXP models, depending on the evolution of 
their equation of state (see Table 1 and Section 2.1). The 
models feature both rapidly and slowly varying equations 
of state and the majority of the models could be classified 
as 'early dark energy' models as they have a non-negligible 
amount of dark energy at early times. 

In order to accurately mimic the dynamics of the origi- 
nal quintessence models at high and low redshift, it is neces- 
sary to use a general prescription for the dark energy equa- 
tion of state which has more parameters than the ubiqui- 
tous 2 variable equation. Parametrisations for w which use 
2 variables are unable to faithfully represent dynamical dark 
energy models over a wide range of redshifts and can lead 
to bi ases when used to constrain parameters (iBassett et al.l 
I2OO4II . Our task has been made easier by the availability 
of parametrizations which accuratel y describe the dynamics 
of the different quintessence models (jCorasaniti fc Copelandl 



Simulations of Quintessential Cold Dark Matter: beyond the cosmological constant 19 



Table 2. The key features in the evolution of the quintessence models simulated. AD{z = 5) is the ratio of the linear growth factor for 
each quintessence model compared to ACDM at 2 = 5. A late time transition in the equation of state is defined as occuring at 2 < 2. 
The AS, CNR, 2EXP and SUGRA models can be considered as 'early dark energy ' models as they have non-negligible amounts of dark 
energy present at early times. 



Model transition type transition rcdshift Qoe{z = 300) AD(z = 5) 



INVl gradual ~ 4.5 ~ 0.009 ~ 50% 

INV2 gradual ~ 5 negligible ~ 10% 

SUGRA rapid ~ 9 ~ 0.01 ~ 20% 

2EXP rapid ~ 4 ~ 0.015 0% 

CNR rapid ~ 5.5 ~ 0.03 0% 

AS rapid ~ 1 ~ 0.11 20% 



I2OO3I : iLinder fc Hutere^bOOSl ). This allows us to modify the 
Friedmann equation in the simulation, using the equation of 
st ate as a function of redshift. W e use the parametrization 
of lCorasaniti fc Copelandl (|2003l ). In its full six parameter 
form, this framework can describe the quintessence model 
back to the epoch of nucleosynthesis. Four parameters are 
sufficient to describe the behaviour of the quintessence field 
over the redshift interval followed by the simulations. With 
this description of the equation of state, our simulations are 
able to accurately describe the impact of the quintessence 
model on the expansion rate of the Universe, from the start- 
ing redshift to the present day. This would not be the case 
with a 2 parameter model for the equation of state. 

In this paper we have taken into account three levels of 
modification from a ACDM cosmology which are necessary if 
we wish to faithfully incorporate the effects of quintessence 
dark energy into a N-body simulation. The first stage is 
to replace the cosmological constant with the quintessence 
model in the Friedmann equation. A quintessence model 
with a different equation of state from w — ~1 will lead to 
a universe with a different expansion history. This in turn 
alters the rate at which perturbations can collapse under 
gravity. The second stage is to allow the change in the expan- 
sion history and perturbations in the quintessence field to 
have an impact on the form of the linear theory power spec- 
trum. The shape of the power spectrum can differ signifi- 
cantly from ACDM on large scales if there is a non-negligible 
amount of dark energy present at early times. This alters the 
shape of the turn-over in the power spectrum compared to 
ACDM. Thirdly, as the quintessence model should be consis- 
tent with observational constraints, the cosmological param- 
eters used for the dark energy model could be different from 
the best fit ACDM parameters. In the three stages of simu- 
lations we look at the effect each of the above modifications 
has on the nonlinear growth of structure. Deconstructing 
the simulations into three stages allows us to isolate specific 
features in the quintessence models which play a key role in 
the growth of dark matter perturbations. 

In the first stage of comparison, in which all that is 
changed is the expansion history of the universe, we found 
that some of the quintessence models showed enhanced 
structure formation at z > compared to ACDM. The 
INVl, INV2, SUGRA and AS models have slower growth 
rates than ACDM. Hence, when normalising to the same erg 
today, structures must form at earlier times in these models 
to overcome the lack of growth at late times. Models such 
as the 2EXP and CNR model have the same recent growth 



rate as ACDM and showed no difference in the growth of 
structure. The difference in linear and nonlinear growth can 
largely be explained by the difference in the growth factor at 
different epochs in the models. At the same growth factor, 
the power in the models only diverges at the 15% level well 
into the nonlinear regime. 

In the second stage, a self-consistent linear theory P{k) 
was used for each quintessence model to generate the initial 
conditions in the simulations. The amount of dark energy 
present at early times will determine the impact on the linear 
theory dark matter power spectrum and the magnitude of 
deviation from the ACDM spectrum. High levels of dark 
energy at early times suppress the growth of the dark matter 
on scales inside the horizon, resulting in a broader turn- 
over in the power spectrum. We found that models with the 
highest levels of dark energy at the last scattering surface, 
such as the AS and CNR models, have linear theory P{k) 
which differ the most from ACDM. The results of the N- 
body simulations of the AS and the SUGRA model show 
a very small increase in nonlinear growth compared to the 
results in Stage I. The increase in the linear theory power 
is on very large scales and does not change the small scale 
growth significantly. 

In our final stage of simulating the effects of 
quintessence, we found the best fitting cosmological param- 
eters for each model, Qmh^, ilhh^ and Ho, consistent with 
current CMB, SN and BAO meaurements. For quintessence 
dark energy models, it is important to consider the changes 
in more than just one cosmological parameter when fit- 
ting to the observational data. For example, for a given 
dark energy equation of state, the values of Qmh^ and Ho 
may change in such a way to compensate one another and 
give similar growth rates and expansion histories to ACDM. 
These compensating effects will be missed if, for example, 
only Slni is changed for the dark energy model as in recent 
work (|Alimi et al.ll2009l ). Models with cosmological param- 
eters which fit the data but were significantly different from 
ACDM were simulated again (Section 14. 3|l . 

We will now summarise and discuss the main results 
for each model. The key features of each of the quintessence 
models are presented in Table 2. The INVl model was un- 
able to fit the data with a reasonable x^/^ (Table A3). This 
toy model had the largest growth factor ratio to ACDM at 
z = 5 and as a result showed the most enhanced growth in 
Stage I of our simulations. The linear growth factor for the 
INV2 model is very different to ACDM at early times and 
gives rise to enhanced growth at z > as seen in Section 
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14.11 This model has negligible dark energy at early times 
and so the spectral shape is not altered in Stage II. In the 
2EXP model the rapid transition to w = — 1 in the equation 
of state early on leaves little impact on the growth of dark 
matter and as a result the power spectra and mass func- 
tion are indistinguishable from ACDM. As both the INV2 
and 2EXP models already agree with cosmological measure- 
ments with very similar values for f2mft^, f2bft^ and Hq to 
ACDM, we did not run these simulations again. 

The SUGRA model has enhanced linear and nonlinear 
growth and halo abundances compared to ACDM at z > 
and an altered linear theory power spectrum shape. The 
mass function results for all stages of our simulations for 
the SUGRA model show enhanced halo abundances &t z > 
0. Analysing the SUGRA power spectra, from a Stage III 
simulation which used the best fit parameters for this model, 
reveals a ~ 5% shift in the position of the first BAO peak. 
We find the distance measure D„ for the SUGRA model 
differs by up to 9% compared to ACDM over the range < 
z < 1.5. Re-scaling the power measured for the SUGRA 
model by the difference in Dv would result in an even larger 
shift in the position of the BAO peaks. 

The CNR model has high levels of dark energy early 
on which alters the spectral shape on such large scales that 
the nonlinear growth of structure is only slightly less than 
ACDM at z < 5. This model has a halo mass abundance at 
z < 5 and BAO peak positions at z = which are the same 
as in ACDM. For z < 0.5 the distance measure, _D„, for the 
CNR model differs from ACDM by ~ 1%, as result there 
would be a corresponding small shift in the BAO peak posi- 
tions. The rapid early transition a.t z = 5.5 in the equation of 
state to KJo = — 1 in this model seems to remove any signal of 
the large amounts of dark energy at early times that might 
be present in the growth of dark matter perturbations. 

The AS model has the highest levels of dark energy at 
early times, and so its linear theory spectrum is altered the 
most. This results in a large increase in large scale power, 
when we normalise the power spectrum to as — 0.8 today. 
The results from Stage III using the best fit parameters show 
both enhanced linear and nonlinear growth at z < 5. The 
linear theory P{k) is altered on scales k ~ O.l/iMpc"^ which 
drives an increase in nonlinear growth on small scales com- 
pared to ACDM. The mass function results in Stage III for 
this model show enhanced halo abundances at z > 0. We 
find that using the best fit cosmological parameters for the 
AS model produces a BAO profile with peak positions simi- 
lar to those in ACDM. At low redshifts there is ~ 1% shift in 
the first peak compared to ACDM after re-scaling the power 
with the difference in the distance measure _D„ between the 
two cosmologies. 

These results from Stage III of our N-body simulations 
show that dynamical dark energy models in which the dark 
energy equation of state makes a late {z < 2) rapid transi- 
tion to Wo = — 1 show enhanced linear and nonlinear growth 
compared to ACDM ai z > and have a greater abundance 
of dark matter haloes compared to ACDM for z > 0. We 
found that dynamical dark energy models can be signifi- 
cantly different from ACDM at late times and still produce 
similar BAO peak positions in the matter power spectrum. 
Models which have a rapid early transition in their dark 
energy equation of state and mimic ACDM after the tran- 
sition, show the same linear and nonlinear growth and halo 
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abundance as ACDM for all redshifts. We have found that 
these models can give rise to BAO peak positions in the 
matter power spectrum which are the same as those in a 
ACDM cosmology. This is true despite these models hav- 
ing non-negligible amounts of dark energy present at early 
times. 

Overall, our analysis shows that the prospects of de- 
tecting dynamical dark energy, which features a late time 
transition, using the halo mass function at 2; > are good, 
provided a good proxy can be found for mass. Parameter de- 
generacies allow some quintessence models to have identical 
BAO peak positions to ACDM and so these measurements 
alone will not be able to rule out some quintessence models. 
Although including the dark energy perturbations has bee n 
found to increase these degeneracies (jWeller fc Lewisll2003l ). 
incorporating them into the N-body code would clearly be 
the next step towards simulating quintessential dark matter 
with a full physical model. Although in many quintessence 
models the dark ene rgy clusters on very l arge scales today 
(fc < 0.02^Mpc"^) (|Weller fc LewisI 120031 ) and the pertur- 
bations are generally small (5de ~ 10~^), these perturba- 
tions may nevertheless have some impact on the dark mat- 
ter structure in a full N-body simulation of the nonlinear 
growth. 
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APPENDIX A: OBSERVATIONAL DISTANCE 
PRIORS 

In this section we outline the method used to find the best fit 
cosmological parameters for each of the quintessence mod- 
els using_CMB^BAO_a^ SN data. The method suggested 
in lKomatsu et akl (|2009l ) employs three distance priors from 
measurements of the C MB together wit h the 'UNION' su- 
pernova samples (jKowalski et al.l l2008t ) and the baryonic 
acoustic oscillations (BAO) in the distribution of galaxies 
(|Percival et al.ll2007l ) to explore the best fit parameters for 
the dynamical dark energy models. In S ect ion [4 . 1 1 and all 
of the quintessence simulations were run using the best fit 
cosmological parameters assuming a ACDM model. While 
this is useful for isolating the effect of the different expan- 
sion histories on the growth of structure, this does not yield 
quintessence models which would automatically satisfy the 
constraints on distance measurements. Using CMB, super- 
novae and BAO data in this way is very useful for test- 
ing and perhaps even ruling out some of the dark energy 
quintessence models. In Section [4.31 we consider the impact 
of using these new cosmological parameters on the nonlinear 
growth of structure. 

These distance priors are derived parameters which de- 
pend on the assumed cosmological model and yield con- 
straints on dark energy parameters which are slightly weaker 
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Table Al. Distance priors based on WMAP observations l lKomatsu et a l for each quintessence model using VL^h? , Q^^h^ and Hq 

parameters from Sanchez ct al. (200s|). These parameters were derived assuming a ACDM cosmology. Z^(2:») is the acoustic scale at the 
epoch of decoupling, z» and R{zt) is the shift parameter. X?otai ~ Xwmap+sn+bao ^^'^ ^ number of degrees of freedom. 





z* Ia{z*) 


R{z*) 




WMAP 5-yr ML 


1090.51 ±0.95 302.10 ± 0.86 


1.710 ± 0.019 





INVl 


261.05 


1.49 


15.34 


INV2 


294.34 


1.67 


1.81 


SUGRA 


284.03 


1.62 


3.88 


2EXP 


303.85 


1.74 


1.09 


AS 


289.69 


1.74 


2.04 


CNR 


306.71 


1.79 


1.37 



Table A2. BA O distance measurements llPercival et al.l^2007^ for each quintessence model using Hm/i^, Q\^h^ and Hq paramet ers from 
ISanchez et al.l (2009). These parameters were derived assuming a ACDM cosmology. A fitting formula proposed by lEisenstein &: Hul 
j 19981 ) was used for the drag redshift ^drag- 





■^drag ^s('2drag) 


rs{zdr.s)/Dv{z = 0.2) 


rsizdr^g)/D^{z = 0.2) 


WMAP 5-yr 


1020.5 ± 1.6 153.3 ± 2.0 Mpc 






Percival et al. f2007) 


154.758 Mpc 


0.198 ± 0.0058 


0.1094 ± 0.0033 


INVl 


152.5 Mpc 


0.208 


0.130 


INV2 


152.7 Mpc 


0.198 


0.121 


SUGRA 


152.5 Mpc 


0.198 


0.121 


2EXP 


152.0 Mpc 


0.192 


0.115 


AS 


143.9 Mpc 


0.183 


0.111 


CNR 


150.7 Mpc 


0.191 


0.114 



than a full Markov Chain Monte Carlo (MCMC) calculation, 
as only part of the full WMAP data is used i.e. the Ci spec- 
trum is condensed into 2 or 3 numbers describing peak po- 
sition and ratios and the polarisation data are ignored. The 
assumed model is a standard FLRW universe with an effec- 
tive number of neutrinos equal to 3.04 and a nearly power 
law primordial power spectrum with negligible primordial 
gravity waves and entropy fluctuations. These WMAP dis- 
tance priors are extremely useful for providing cosmological 
parameter constraints at a reduced computational cost com- 
pared to a full MCMC calculation. 

We shall briefly review the distance scales used in this 
paper and the method for finding the best fit parameters for 
the dark energy models. From measurements of the peaks 
and troughs of the acoustic oscillations in the photon-baryon 
plasma in the CMB it is p ossible to measure two distance 
ratios (iKomatsu et al.ll2009l ). The first ratio is quantified by 
the 'acoustic scale', Ia, which is defined in terms of the sound 
horizon at decoupling, rs{z,) and the angular diameter dis- 
tance to the last scattering surface, Da{z,), as 



Ia = (l + z,) 



■kDa{z,) 



(Al) 



Assuming a flat universe, the proper angular diameter dis- 
tance is defined as 



"""^^^ =(1^)70 HiP-y 

and the comoving sound horizon is given by 



r»(z) 



73 



1/(1+^) 



da 



a2J/(a)^l + (3nb/4n^ 



(A2) 



(A3) 



where O-y = 2.469 x lO'^h-^ for Tcmb = 2.725K 
(jKomatsu et al.| |2009l ) and ilb is the ratio of the baryon en- 
ergy density to the c ritical density. We shal l use the fitting 
formula proposed bv lHu fc Sugivamal l|l996l ) for the decou- 
pling epoch Zt which is a function of Q,hh^ and flmh^ only. 
The second distance ratio measured by the CMB is called 
the 'shift parameter' (|Bond et al.lll997l ). This is the ratio of 
the angular diameter distance and the Hubble horizon size 
at the decoupling epoch which is written as 



-z,)Da{z, 



(A4) 



Eg. I A4I assumes a standard radiation and matter dominated 
epoch when calculating the sound horizon. The expression 
for the shift parameter will be modified for quintessence 
models of dark energy. The proper expressi on for the shift 
parameter is given by (iKowalski et al. 



cpressio 

lliooi) 



J?(z,) = i?std(z») 



d 



:/V»n.(l + 2)3 

SZ dzHo/H{z) 



(A5) 



where is the standard shift parameter given in Eq. IA4I 
This correction to the shift parameter can be substantial 
for quintessence models with non-negligible amounts of dark 
energy at early times and so we include this correction for all 
of the scalar field models in this paper. The 5-year WMAP 
constraints otiIa, R and the redshift at decoupling are the 
WMAP distance prio rs used to test models of dark energy 
jKomatsu et al.ll2009l ). 

The angular diameter distance at the decoupling epoch 
can be determined from measurements of the acoustic oscil- 
lations in the CMB. These baryonic acoustic oscillations are 
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Table A3. Best fit values for Q^h-^, Oj^/i^ and Hg with 68.3% confidence intervals from minimising Xtotai ~ XwMAP+sci±aAQ.j2I_22£il 
quint essence model. wCDM WMAP 5-year are the parameter constraints assuming a dynamical dark energy model jKomatsu et al.l 
I2OO9I) . 

102Qb/i2 Ho (km/s/Mpc) Qn-^h^ xL.J'' 



ACDM WMAP 5-yr Mean 2.267 70.5 ±1.3 0.1358 lo.oo36 

wCDM WMAP 5-yr Mean 2.27 ±0.06 69.7 ±1.4 0.1351 ±0.0051 

Sanchez et al. (20091 2.267 tjjlos'^ 71.5 ±1.1 0.13343 ±0.0026 1.09 



INVl 


3.78 ±0.145 


63.13 ±0.5 


0.115 ±0.0103 


2.27 


INV2 


2.35 ±0.094 


68.21 ±0.7 


0.124 ±0.0065 


1.07 


SUGRA 


2.68 ±0.105 


67.63 ±0.7 


0.111 ±0.0075 


1.25 


2EXP 


2.22 ±0.115 


70.01 ±0.8 


0.138 ±0.0031 


1.05 


AS 


2.12 ±0.121 


70.42 ±0.9 


0.086 ±0.0121 


1.07 


CNR 


2.09 ±0.185 


70.05 ±1.2 


0.140 ±0.0133 


1.12 



Table A4. WMAP distance priors jKomatsu et al.ll2009[ ) for each quintessence model using the best fit parameters Clmh'^, 0^/1^ and 
Ho given in Table A3. 





z* 


Ia{z*) 


R{z*) 


ACDM WMAP 5-yr ML 


1090.51 ±0.95 


302.10 ± 0.86 


1.710 ± 0.019 


Sanchez et al.2009 


1090.12 ± 0.93 


301.58 ± 0.67 


1.701 ± 0.018 


INVl 


1076.17 


292.54 


1.519 


INV2 


1088.71 


301.69 


1.676 


SUGRA 


1083.96 


298.51 


1.596 


2EXP 


1091.75 


302.91 


1.749 


AS 


1087.98 


300.23 


1.684 


CNR 


1093.97 


303.51 


1.809 



also imprinted on the distribution of matter. Using galaxies 
as tracers for the underlying matter distribution the cluster- 
ing perpendicular to the line of sight gives a measurement of 
the angular diameter distance, Da{z)- BAO data also allow 
us to measure the expansion rate of the universe, H{z), from 
observations of c l usteri ng along the line of sight. Recently, 
iGaztanaga et all l|2008l ) made a direct measurement of the 
Hubble parameter as a function of redshift providing for the 
first time a measure of Da{z) and H{z) individually. Us- 
ing a spherically averaged correlation function to reveal the 
BAO signal results in a n effective distance measure given by 
lEisenstein et aLlbOOsI ) 

Dv{z) = (^(1 + zfDliz)-^^ . (A6) 

It is the ratio of Dv{z) to the sound horizon, r^, at the 
drag epoch, Zdrag, which determines the peak positions of 
the BAO signal. The drag epoch is the redshift at which 
baryons are separated from photons and is slightly later 
than the decoupling epoch, z*. For a wide angle survey, 
Dv is used, which is motivated on dimensional grounds and 
equal sampling of all axes (e.g. for a penc il beam survey 
would have different exponents of Da and J/V lPercival et ahl 
((2003) provide rs{zd)/Dv{z) at two redshifts, z = 0.2 and 
z = 0.35, taken from the Sloan Digital Sky Survey (SDSS) 
and Two Degree Field Galaxy Redshift Survey (2dFGRS). 
The two values are rs{zd) / Dv{G.2) = 0.198 ± 0.0058 and 
rs(zd)/Dv{Q-ib) = 0.1094 ± 0.0033. 



The UNION supernovae compilation (|Kowalski et al.l 

l2008h consists of 307 low r edshift SN all p rocessed using 
the SALT Ught curve fitter ||Guv et al.1 120051 '). This compil- 
ation includes older data sets from the Supernova Legacy 
Survey and ESSENCE Survey as well as a recent dataset 
observed with HST. Type la supernovae data is extremely 
useful in breaking parameter degeneracies such as the w, 
SIde degeneracy in the CMB data. A wide range of these 
two parameters can produce similar angular diameter dis- 
tances at the redshift of decoupling and so SN constraints, 
which are almost orthogonal to CMB constraints, help to 
reduce this parameter space. The current SN data cover a 
wide range of redshift, 0.02 55 2: ^ 1.7, but is only able to 
weakly constrain a dynamical dark energy equation of state, 
w, ai z ^ 1. Also, due to a degeneracy with Qui, the current 
SN data by themselves are not able to tightly constrain the 
present value of w and including measurements involving flm 
such as CMB or BAO observations break this degeneracy. 

Following the prescription of iKomatsu et al] (|2009l ) for 
using the WMAP distance priors it is necessary to find 
the vector x = {lA,R,Zt) for each quintessence model in 
order to compute the likelihood, C, as — —2\njC — 
(Xi — dijC^j (Xj — dj), where d = {La ,2, ) 

and C~:^ is the inverse covariance matrix for the WMAP 
distance priors. 

In order to find the best fit cosmological parameters for 
each quintessence model we minimise the function Xtotai = 
XwMAP + xIao + XsN with respect to Q.mh^ , Q.hh' and Hq. 
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Table A5. BAO distance measurements dPercival et al.l^2007^ for each quintessence model using the best fit parameters Vl^h^ , Q.^h'^ 
and Hq given in Table A3. 





^drag 


r-sCZdrag) 


rs/Dv{z = 0.2) 


rs/Dv{,z = 0.35) 


WMAP 5-yr 


1020.5 ± 1.6 


153.3 ± 2.0 Mpc 






Percival et al. f2007) 




154.758 Mpc 


0.198 ± 0.0058 


0.1094 ± 0.0033 


INVl 


1045.1 


146.259 Mpc 


0.1765 


0.1103 


INV2 


1021.2 


154.946 Mpc 


0.1921 


0.1167 


SUGRA 


1026.4 


155.803 Mpc 


0.1908 


0.1161 


2EXP 


1019.9 


150.983 Mpc 


0.1879 


0.1123 


AS 


1010.5 


157.745 Mpc 


0.1947 


0.1161 


CNR 


1017.1 


150.597 Mpc 


0.1876 


0.1128 



In appendix D of iKomatsu et al.l (|2009l ) it can be seen that 
including the systematic errors has a very small effect on 
the ACDM parameters but can have a significant impact on 
dark energy parameters. Usi ng a two parameter e quation 
of state for the dark energy iKomatsu et al.l l)2009l ) found 
that the parameter constraints weakened considerably after 
including systematic errors. In calculating Xsn in this paper 
we have used the covariance matrix for the errors on the SN 
distance moduli without systematic errors. 

Table Al shows the WMAP distance priors computed 
for each d ark energy model usi ng the cosmological param- 
eters from ISanchez et al.l (|2009l ) . The BAO scale and drag 
redshift, Zd, are given in Table A2 using the same param- 
eters. From these tables it is clear that some quintessence 
models with ACDM cosmological parameters fail to agree 
with the distance measurements within the current con- 
straints. 

With the assumption that f2mft^, 0.hh^ and Ho are 
tightly constrained by WMAP, BAO and SN data, and as a 
result their posterior distribution is close to a normal distri- 
bution, minimising x?otai = Xwhap+xIao+xIn with respect 
to these three parameters will be the same as marginalising 
the posterior distribution. We have fixed the dark energy 
equation of state parameters for each quintessence model 
and the 68.3% confidence intervals for each parameter from 
minimising Xtotai are shown in Table A3. The final column 
in this table is /v where v is the number of degrees of 
freedom. From Table A3 it is clear that the INVl model is 
unable to fit the data and has a poor /v = 2.27 statis- 
tic. Most of the quintessence models favour a lower Q.^h?' 
compared to ACDM in order to fit the distance data. As 
can be seen from Table A3 the confidence intervals on the 
three fitted parameters Q,mh? , f2b/i^ and Hq are quite large. 
Once the best fit parameters from Table A3 are used, all 
of the quintessence models apart from INVl which we rule 
out, produce a better fit to the data, as seen in Tables A4 
and A5, for the WMAP distance priors and the BAO dis- 
tance measures respectively. As we noted earlier the WMAP 
distance priors do not contain all of the WMAP power spec- 
trum data and only use the information from the oscillations 
present at small angular scale (high multipole moments). Ne- 
glecting the Sachs- Wolfe (SW) effect at large angular scales 
(small multipole moments) as well as polarisation data lead 
to weaker constraints on cosmological parameters in these 
dark energy models. We have not considered how these dis- 
tance priors wo uld change wit h the inclusion of dark energy 
perturbations (|Li et al.l l2008l ). These results are in agree- 



ment with previous work fitting cosmological parameters of 
quintesse nce models using WMA P first year CMB data and 
SN data (jCorasaniti et al.ll2004l ). 

This paper has been typeset from a TeX/ ffl^jX file prepared 
by the author. 



